Gapped state in three-leg S — I Heisenberg tube 
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We study the ground state of three-leg S = i Heisenberg tube using the density-matrix renormal- 
ization group method. The dimerization order-parameter and spin-excitation gap are calculated in a 
wide range of leg exchange interactions. We confirm that a gapped state with the dimerization order 
is realized even when the leg exchange interactions are ferromagnetic. Furthermore, the topological 
configuration of spin-singlet pairs in the ordered state is determined from the results of the Berry 
phase of each spin coupling as well as the structure factor of singlet-singlet correlation functions. 
We find that there exist three kinds of configurations of the valence-bond states depending on the 
ratio of leg and rung exchange interactions. 
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Odd-leg spin ladder belongs to the same universality 
class as single chain; thus, the ground state is compre- 
hended as a gapless spin-liquid (or a Tomonaga-Luttinger 
liquid) X However, if the periodic boundary conditions are 
applied in the rung direction, i.e., a tube is shaped, the 
spin states are dramatically changed. It is due to the oc- 
currence of a geometric property called frustration, which 
is today a hot topic in condensed matter physics.— In gen- 
eral, it would appear that the fundamental low-energy 
physics of any odd-leg spin tube is essentially epitomized 
by that of three-leg one. So far, it has been recognized 
that, as long as all the exchange interactions are antifcr- 
romagnetic, the three-leg spin tube can be spontaneously 
dimerized to avoid (or to reduce) the frustration and the 
spin excitations arc gapped£i^ 

An ideal nanotubular material with odd number of legs 
is vanadium oxide Na2Vs07, which may be regarded as 
a S = h nine-leg Heisenberg spin tube system^ In ex- 
periments (Ref. 0), the 23 Na NMR response, the dc- and 
ac-magnetic susceptibilities, and the specific heat reveal 
that above 100 K the system is considered as param- 
agnetic; whereas, below 100 K most of the localized V 
magnetic moments (S — \) form a collection of spin- 
singlet dimers with gaps A ~ — 350 K and the re- 
maining small fraction of them forms spin-triplet bound 
states with gaps A ~ — 15 K; and, the degeneracy of 
the triplet ground states is lifted by a phase transition 
at 0.086 K. The mechanism of the gap opening is still 
open issue. Besides, this material has attracted consid- 
erable attention from a standpoint of entanglement in 
spin systems Moreover, it is of great interest to seek 
a relevance to a spin-liquid state observed in a three-leg 
5=| spin tube system CsCrF^ 

The low-energy spin Hamiltonian of Na2 V3O7 has been 
proposed by some theoretical groups, nevertheless, it is 
still controversial. Both the rung (J±) and leg (J||) ex- 
change interactions seem to be very sensitive because 
the estimated values are quite different, even in ferro- 
and antiferromagnetic characteristics, depending on the 



approaches. An overview of the results are as follows, 
(i) The ab initio microscopic analysis^ both J± and Jy 
are antiferromagnetic, but Ju is frustrated and the mag- 
nitude is much smaller than J±; (ii) the first-principle 
calculations:— J± is antiferromagnetic, while Jy is ferro- 
magnetic, and they have the same order of magnitude; 
and, (hi) the first-principles density functional theory^ 
J|i is ferromagnetic and Jj_ is ferro- or antiferromagnetic. 
Possibly, the point to be grasped next is whether the 
dimerization order with finite spin gap can occur when 
ferromagnetic exchange interactions are contained in the 
odd-leg spin tube. 

We thus consider the ground state and low-lying ex- 
cited states of three-leg S — \ Heisenberg spin tube and 
provide new insights especially for the case that the fer- 
romagnetic exchange interactions are taken into account. 
The Hamiltonian is given by 

3 L L 
H = J|| S a j ■ S a ^ + i + J± S a ^ ■ S a '^,(l) 

a=l i=l a(=jta') i=l 

where S a ^ is a spin-i operator at rung i and leg a. We 
here restrict the rung interactions J± to be antiferromag- 
netic because it is obvious that no dimerization order 
occurs with ferromagnetic rung interactions. The leg in- 
teractions J|| are varied from ferro- to antiferromagnetic 
range. Although several theoretical studies have been 
carried out on this and similar models, only antiferro- 
magnetic interactions are considered*^— . 

In this paper, we calculate the dimerization order- 
parameter and the spin-excitation gap to clarify in which 
range of the exchange interactions the ordered state is re- 
alized. Also, the Berry phase of each coupling and the 
structure factor of singlet-singlet correlation functions 
are calculated to check topological configuration of the 
spin-singlet pairs in the ordered state. For those calcula- 
tions, the density- matrix renormalization group (DMRG) 
technique^ is applied. We investigate tubes with several 
kinds of length up to L = 312, i.e., 312 x 3 cluster, under 
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FIG. 1: (a) Dimerization order-parameter D extrapolated to 
the thermodynamic limit £^ooasa function of J±/J«. (b) 
Extended figure of (a) for —3 < J±/J\\ < 1. 



the open boundary conditions (OBC) in the leg direction, 
unless otherwise stated. The density-matrix eigenstates 
up to m = 2400 are kept in the renormalization pro- 
cedure and all quantities are extrapolated to the limit 
m — > oo. In this way, the discarded weight is less than 
1 x 10~ 7 , while the maximum error in the ground-state 
energy is less than 10 — 10~ 6 . 

Let us first evaluate the dimerization order-parameter 
in a wide range of Jj_/J|| for exploring the presence or 
absence of long-range dimerized state. The order pa- 
rameter is featured by an alternation of nearest-neighbor 
spin-spin correlations, S(i) — — (o a ^-S a ^+i), where (■ • ■) 
denotes the ground-state expectation value. As the OBC 
breaks translational symmetry in our calculation, the 
dimerized state is observed as a Friedel oscillation. Gen- 
erally, the amplitude of the Friedel oscillation at system 
center decays as a function of system size. In case that 
it persists for arbitrarily large system size, we can judge 
that a long-range order exists. Thus, the dimerization 
order-parameter is defined as 



D= lim \S(L/2)-S(L/2 + l) 



(2) 



Nonzero value of D indicates the presence of long-range 
dimerization order in our model. In Fig. [T] (a), the re- 
sults of D are shown as a function of J±/ J\\ ■ We find that 
the dimerization order appears over the entire region of 
J±/J\\(¥ 1 0)- For both positive and negative J±/J\\ val- 
ues, D increases rapidly at J±/\J^\ ^ 10 and keeps al- 
most constant at J± > 10: in the limit of Jj_/J\\ — > oo 
(— oo), it is saturated to D ~ 0.0673 (0.0183). However, 
quite different behaviors are seen if we take a closer look 
at small J_l/\J\\\ regime [see Fig. ffl(b)]. At Jj_/J\\ > 0, 
D is discontinuously enhanced when J± is switched on, 
then goes through a minimum around J± = 0.1, and 
increases almost linearly from Jj_ « 0.2 to 5; while at 



J±/J\\ < 0, D increases gradually with increasing J±, 
like D ~ exp[-0(l/J]_)]. 

The dimerization order implies a complete occupation 
of the system with spin-singlct pairs. Which poses a ques- 
tion on the topological configuration of the pairs. In or- 
der to solve it, we investigate the quantized Berry phase 
in the ordered state. The Berry phase is defined by 
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I I A(<j>)d<t>, 





(3) 



where A{<j>) is the Abelian Berry connection, A{4>) = 
(Vvl^VvO w ith the ground state The Berry phase 

is quantized as or tt (mod 27r) if the system has spin gap 
during the adiabatic continuation and time reversal sym- 
metry; and "undefined" if a gapless excitation exists. We 
introduce a local perturbation by a twist of the nearest- 
neighbor connection, § a>i ■ S a ,j -> |(e _l *5^ i 5~, j + 

ei *S-iS+, tj ) + SZ ti SZ,j. We here pick up a dimerized 
pair of triangles, i.e., including six spins, and evaluate the 
Berry phases of the leg bond (Ticg) for a = a' ,j = i + 1 
and of the rung bond (7 rU ng) for a ^ a',j = i. Note 
that the dimerized pair of triangles must include three 
spin-singlet pairs. 

In Fig. [2 the Jj_/ Jy-depende nee of 7icg and 7rung 
are shown and the corresponding configurations are also 
schematically described. We find that there are three 
kinds of the spin-singlet configurations along Ji/Jy. Let 
us now see the case that J|| is antiferromagnetic. A recon- 
struction of the valence bonds is seen at J±/J\\ ~ 5. This 
would correspond to the crossover between constant-A 
and proportional- A regions around Jx/J\\ = 5^1 For 
J±/J\\ ^ 5, the Berry phases of leg and rung bonds are 
denoted by tt and 0, respectively. It means that all singlet 
pairs are formed on the leg bond [Fig. [2(iii)]. Whereas, 
for < J±/J\\ 5, both 7i cg and 7 run g are denoted by 
0. It may be interpreted if we assume that both the leg 
and rung bonds are involved to form spin-singlet pairs: 
one pair is formed in cither one of three legs and the 
other four spins form a couple of pairs in the two rungs 
[Fig.[2][ii)]. This valence-bond state is threefold degener- 
ate, so that the spin-singlct pairs would be not detected 
as local ones. Then, we turn to the case that Jy is ferro- 
magnetic. A valence-bond state is detected on the rung 
bond and a gapless excitation is found in the leg bond; 
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FIG. 2: Berry phases on the rung bond (7 run g) and leg bond 
(Ticg) along the parameter J±/Ju. Schematic pictures of the 
valence-bond state are also shown. Two dots linked by bold 
line denote a spin-singlet pair. 
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FIG. 3: (a) Spin gap A extrapolated to the thermodynamic 
limit L — > oo as a function of J±/J«. (b) Extended figure of 
(a) for -3 < Jl/J|| < 1. 




however, the spin excitations are gapped in bulk limit as 
shown below. It will be settled if we suppose that two 
spin-singlets are formed on rungs of each triangle and 
the other one is between two sites on different legs and 
triangles [Fig. |2{i)]. We call the last pair "diagonal spin- 
singlet" and the effective antiferromagnetic exchange in- 
teraction is denoted by Jd hereafter. 

And now, of particular interest is the magnitude of 
spin-excitation gap especially in the ferromagnetic Jii 
case. It will be clarified by calculating an energy differ- 
ence between the singlet ground state and the first triplet 
excited state, 

A = lim [E(L, 1) - E(L, 0)], (4) 

where E(L, S z ) is the ground-state energy of the system 
of length L with the z-component of total spin S z . We 
note that the system length must be taken as L = 21, 
with Z(> 1) being an integer to maintain total spin of 
the ground state as S = 0. In Fig. [31 we show the re- 
sults of A as a function of Jx/Jn. We see that the gap 
opens except at Jj_ = 0, as anticipated from the results 
of dimerization order-parameter. For both J|| > and 
Jii < 0, roughly speaking, A starts to increase propor- 
tionally to Jx and shift into almost constant for larger 
Jx- It means that the lowest singlet-triplet excitations 
for small and large J_l/|J|i| are scaled by distinct ex- 
change interactions. When Jm is antiferromagnetic, as 
suggested in our previous paper^i A is approximately 
scaled by Jx (J\\) in the small (large) Jx/J\\ regime. On 
the other hand, A seems to be always scaled by Jd since 
the diagonal spin-singlet pair has the smallest binding 
energy when J|| < 0. It may be interpreted as follows: 
in the perturbative sense, Jd is proportional to Jx (| J\\\) 
for Jx J|| {Jx J||); accordingly, A is scaled by 
Jx (|J|||) in the small (large) Jj_/|Jii| regime as in the 
case of antiferromagnetic Jy. In fact, A for J±/J\\ < 



FIG. 4: Singlet-singlet correlation functions 5?(k) for (a) 
Jy =1,J ± = 1, (b) Jy = 1, J x = 10, and (c) J y = -1, J± = 
10. 

is about a factor of 6 smaller than that for Jx/J\\ > 0. 
As a consequence, we obtain A = 0.0437 in the limit of 
Jj_/J|| = -oo; A = 0.254 in the limit of Jj_/J|| = oo. We 
note that the Jj_-dependence of A looks similar to that 
of D but except when Jx/\ J\\ \ is very small. 

Lastly, in order to determine the periodicity of align- 
ment of the valence-bond configurations, we calculate 
singlet-singlet correlation functions 

= JlY. ex P[^ ' (Vi - r aljj )\ (5) 

aa'ij 

with spin-singlet number operator 

'jL.i — ^ $a,i ' $a',ji (6) 

where (a',j) = (a, i + 1), (a + 1, i), and (a + 1, i + 1) 
are chosen for the leg- leg (LL), rung-rung (RR), and 
diagonal-diagonal (DD) correlations, respectively. Here, 
the system size is fixed at L = 24 using the periodic 
boundary conditions. Thus, relations a + 3 = a and 
i + 24 = i arc fulfilled. 

In Fig. HI the correlation functions ^(k) at (a) Jii = 
1, Jx = 1, (b) Jy = 1, J± = 10, and (c) J\\ = -1, Jj_ = 10 
for L = 24 are shown [the deviation from those for 
L = 12 has been confirmed to be negligibly-small]. We 
initially point out that large values of ^(tt, |7r) in the 
LL and of =5^(0, |7r) in the RR correlations jointly in- 
dicate a straight link of the state (ii) in Fig. [21 which 
is described as alignment (II) in Fig. [5J whereas, large 
value of ^(ir, 0) in the LL correlation suggests a link 
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FIG. 5: Possible configurations of spin-singlet pairs for re- 
gions (I) J ± /J|| < 0, (II) < Jx/Jy < 5, and (III) 5 < Jx/Jy. 
Two dots linked by bold line denote a spin-singlet pair. 

of state (iii), which corresponds to alignment (III). For 
J|| = 1, J± = 1, therefore, it would appear that the align- 
ments (II) and (III) coexist [see FigQJa)]; however, the 
alignment (III) must be less dominant considering the 
results of the Berry phase. As J± increases, in the LL 
correlation S*(tt, |7t) is rapidly diminished and ^(ir, 0) is 
rather enhanced [see Figj4jb)]; as a result, the alignment 



(III) becomes dominant for large J±/J\\ region. Then, 
we turn to the case of J|| < [see FigQJc)]. As expected, 
^(k) in the LL correlation is always zero indicating no 
spin-singlet pair in the leg bond; instead, the DD corre- 
lations are significant. The enhancements of =5^(7r, |7t) in 
the RR and S^{tt^tt) in the DD correlations are seen. It 
indicates a straight alignment of the state (i) , as shown in 
Fig. [5] (I). We argue that a zigzag chain, which is denoted 
as a dotted line in Fig. [SJ is dimerized to remove the frus- 
tration instead of being dimerized in the leg direction for 
antiferromagnetic J\\ . 

In conclusion, we study the ground-state properties of 
three-leg S = ^ Heisenberg tube with antiferromagnetic 
rung exchange interactions. Using the DMRG method, 
the dimerization order-parameter, spin-excitation gap, 
Berry phase, and structure factor of singlet-singlet cor- 
relation functions are calculated. We confirm that the 
dimerization order occurs for both ferro- and antiferro- 
magnetic leg exchange interactions, and the spin excita- 
tions are always gapped. Also, we find that there are 
three kinds of phases according to topological configura- 
tion of the spin-singlet pairs: (I) diagonal-singlet regime 
for J±/J\\ < 0, (II) rung-singlet regime for < J±/ < 
5, and (III) leg-singlet regime for J±/Jn > 5. 
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